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We investigate the Gamow- Teller strength distributions in the electron-capture direction in nuclei 
having mass A = 90 — 97, assuming a ®**Sr core and using a reahstic interaction that reasonably 
reproduces nuclear excitation spectra for a wide range of nuclei in the region as well as experimental 
data on Gamow- Teller strength distributions. We discuss the systematics of the distributions and 
their centroids. We also predict the strength distributions for several nuclei involving stable isotopes 
that should be experimentally accessible for one-particle exchange reactions in the near future. 
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I. INTRODUCTION 

New frontiers of nuclear structure experiments to 
probe the Gamow- Teller distributions in medium-mass 
nuclei are currently being pursued. These experiments 
will be able to measure Gamow- Teller data in the mass 
90-100 region. Extensive theoretical studies have been 
devoted to Gamow- Teller total strengths and strength 
distributions in Is-Od shell nuclei (mass A = 16-40 nu- 
clei) [1] and the 0/-lp shell (mass A = 40-80 nuclei) [2-4]. 
Due to an excellent agreement between shell model re- 
sults and the available experimental data, the calculated 
results have been used extensively to predict numerous 
Gamow- Teller strength distributions in nuclei that have 
not yet become experimentally accessible [5]. 

In addition to their nuclear structure interests, an ap- 
propriate description of Gamow- Teller transitions in nu- 
clei directly affects the early phases of type II supernova 
core collapse since electron capture rates are partly de- 
termined by them. The effects of the improved rate es- 
timates are rather dramatic, as was recently discussed in 
Refs. [6, 7]. In addition to the standard Gamow- Teller 
transitions, first- and second-forbidden transitions con- 
tribute to the electron capture rates in the supernova 
environment. For terrestrial experiments, the primary 
focus is on the Gamow- Teller transitions. 

Recently, Zegers et al [8] proposed to measure the 
Gamow- Teller distributions using stable Zr and Mo iso- 
topes as targets in (t, ^He) reactions [9]. Estimates 
indicate that the Gamow- Teller strength is sufficiently 
large to be measured. In this paper, we will investi- 
gate these transitions using standard shell-model diago- 
nalization techniques for 36 nuclei with the mass number 
90 < A < 97 (Z = 40-47, N = 50-57). To validate the in- 
teraction, we also studied excitation spectra in those and 
other nuclei in the region. Since our model space does 
not contain all spin-orbit partners, i.e. it is not a com- 
plete Qhio calculation, the total Gamow- Teller strength 
will be overestimated in our calculations. We adopt a 
single quenching factor similar to the one discussed in 
Ref. [10]. We estimated this factor based on recent ex- 
perimental data on ^^Ag [11]. We used this measurement 
to gauge our calculation for two reasons. First, it used 



the total absorption spectrometry which accounts also 
for the weak 7-ray cascades that follow the P'^ decay. 
Second, almost all total Gamow- Teller strength is inside 
the Q-window. We note that this factor need not be 
universal as it is simply a phcnomcnological tool at this 
point. 

The remainder of this paper is organized as follows. 
In Section II, we present results on the nuclear spectra, 
generated with an effective interaction that uses ^®Sr as 
a core, and compare them to experiment. In section III, 
we present our shell-model diagonalization results for the 
Gamow- Teller strength distributions and compare to ex- 
periment when available. We also present systematics of 
the Gamow- Teller centroids. In section IV, we discuss 
the distributed-memory shell-model computer code that 
we developed and used for these calculations. Finally, we 
conclude and give a perspective in Section V. 



II. 



CALCULATED SPECTRA USING THE ^^SR 
CORE 



We perform our shell-model diagonalization calcula- 
tions in a model space taking *^Sr as the core nucleus 
and allowing excitations within the valence space of Ipi /2 
and 0^9/2 proton shells and 1^5/27 2si/2, 1^3/2, 037/2, and 
O/111/2 neutron shells. While it cannot be used for calcu- 
lations of /3-decays, it appears suitable for Gamow- Teller 
distributions in the electron-capture direction. The ef- 
fective interaction [12] was derived from a CD-Bonn po- 
tential [13] using the machinery of many-body pertur- 
bation theory [14]. We use the following single-particle 
energies: £(^1/2) = 0.0 MeV and £(39/2) = 0.9 MeV 
for protons; and £(^5/2) = 0.0 MeV, £(51/2) — 1-26 
MeV, £(^3/2) = 2.23 MeV, £(57/2) = 2.90 MeV, and 
^(^11/2) = 3.50 MeV for neutrons. A slightly different 
version of this interaction was used to describe Sr and 
Zr isotopes [15]. We have not attempted to adjust the 
interaction to obtain a better fit to experimental data 
[16]. 

We calculated low-energy spectra of more than 50 nu- 
clei with masses 90 < A < 98, 38 < Z < 48, and 
50 < A^ < 58. General agreement between the calcu- 
lated lowest states and experimentally observed states 
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FIG. 1: Experimental [17] and calculated spectra in '^''Pd, ^''Rh, ^''Ru, ^''Tc, ^"Mo, and ^''Nb. Experimental levels are shown in 
the left column; the middle and right columns display calculated negative and positive parity states, respectively. Experimental 
positive parity states are shown using dashed lines; negative parity states are shown using dotted lines; states with known spin 
and uncertain parity are shown using dash-dotted lines; the states having undetermined both spin and parity are shown using 
a shorter dashed line. The sign '&' indicates that assignments refer to different experimental states, with the second one lying 
at a higher excitation energy. Since there is no spin ambiguity in the calculation, a comma is used in the middle and right 
columns in a similar situation. In some cases a subscript number is used to distinguish levels having the same spin and parity. 
At most, 10 calculated negative parity states are shown for each nucleus. Only 1+ calculated states are shown in ^^Rh above 
Ex « 1.1 MeV. Only negative parity states and positive parity states with J > 5 are shown for ®®Ru above Ex = 2.5 MeV in 
experimental spectrum. Experimental positive parity states are cut-off ^^Mo after Ex = 2.5 MeV. 
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FIG. 2: Experimental and calculated spectra in ®®Zr, ®^Y, and 
^®Sr. For notations, see caption of Fig. 1. Excitation energy 
of the isomeric state in ®^ Y with J^ = 8+ is not known [17] . 



is satisfactory. We judged the agreement based on re- 
production of low-lying states up to a chosen excitation 
energy. For odd-odd isotopes, which have higher density 
of states, the upper limit was chosen to be 1 MeV. For 
even-even nuclei the limit was up to 3 MeV. Not all ob- 
served states were found in the model space, as would 



be expected from a restricted calculation, and for some 
nuclei our calculations suggested some low-lying states 
that have not yet been observed. 

The interaction generally reproduces the correct spin 
for the lowest states of both parities as well as their order- 
ing, though there are cases where some levels are inter- 
changed. The energy splitting between the lowest states 
with different parities is reproduced with varying success, 
although this is difficult to judge for some nuclei because 
of the lack of experimental information. 

For representative spectra which indicate the overall 
quality of the interaction, we show nuclei having mass 
A = 96 in Figs. 1-2. The maximum energy range shown 
in the plots is varied following the density of states. From 
these figures we observe that the spectra of even-even nu- 
clei (96pd, 96Ru^ 96mo, s^Zr, and ^^Sr) are reproduced 
well. For ^^Pd there are more calculated states than ex- 
perimentally known. In some nuclei the model space is 
insuflRcient to describe all observed states. Odd-odd nu- 
clei (^^Rh,^^Tc, ^^Nb, and ^®Y), having more states, are 
also more difficult to describe although, even here, the 
interaction performs reasonably well. The position of 8"*" 
isomer in ^^Y is not known experimentally [17]. This 
state appears in the calculation at a relatively high exci- 
tation energy, 1.1 MeV above the lowest positive parity 
state which was calculated to be 5+. The nucleus ^^Tc 
reflects a situation where the lowest states are experimen- 
tally very close (in this nucleus there are 6 states in the 
energy range of 50 keV), while the calculation reproduces 
the states but not their energies (the calculated range is 
310 keV). A similar situation occurs in ^^'^^Nb. 

These ^ = 96 nuclei reflect the situation in other cases 
as well, with a general conclusion that the interaction re- 
produces excitation spectra reasonably well, though fine 
tuning might increase the accuracy. We do not discuss 
them in a greater detail, since the focus of our paper is 
Gamow- Teller distributions. 



III. GAMOW-TELLER STRENGTH 
DISTRIBUTIONS 

Our study focuses on the Gamow- Teller transitions 
from the lowest positive parity states which is natural for 
the most nuclei in the region above *^Sr, with the excep- 
tion of the Y isotopes where the odd proton in the Pii2 
shell is responsible for low- lying negative parity states. 
Since our model space is not sufficient to reasonably re- 
produce negative parity states in Sr isotopes where no 
valence protons are available, we do not calculate the 
transitions between these two isotope chains. Among the 
calculated nuclei, there are three cases where we chose 
the lowest experimental state to be the initial state for 
GT excitations rather than using our calculated lowest 
energy state. This affected two A^ = 51 nuclei, ^^Nb and 
^^Tc, where the calculation places 2+ to be the lowest 
state, and the nucleus ^^Tc. 

The Gamow- Teller strength was calculated using the 



formula 



GT+ = {cTTf = 
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(1) 

To obtain the strength distribution, we used the method 
of moments [18]. We performed 33 iterations for each 
Jf in all nuclei except for the decays of ^^Mo, where 
we did 24 iterations per final state, and ^^Ag, where a 
complete convergence was achieved. The GT+ strength 
inside the experimental Qec window [19] is marked as 
^Qec- This value is only an estimate, since we did not 
strive to achieve the convergence of states inside the Q- 
window. 

As discussed above, our calculation overpredicts the 
Gamow- Teller strength; thus we include a hindrance fac- 
tor, h, so that S{GT+) = GT+/h. This factor is found by 
comparing experimental data to the calculated Gamow- 
Teller total strength. For nuclei around ^°"Sn, the single- 
particle estimate of the Gamow- Teller strength is com- 
monly used, since the main contribution comes from a 
transition of a 59/2 proton into a gi/2 neutron. The esti- 
mate is given by a formula (see e.g. [20]): 
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where iVg/2 is the occupation of the 179/2 shell by protons, 
and N'Y/2 is the occupation of 37/2 shell by neutrons in 
the initial state of a parent nucleus, and GT+(^'^'^Sn) = 
17.78 is the single particle estimate of the total Gamow- 
Teller strength for ^"^Sn. In the simplest, non-interacting 
shell model, the occupation numbers are replaced by the 
numbers of valence particles in the corresponding shells. 
This simplest estimate does not exactly reproduce our 
calculated strength even though the values are close. It 
underestimates the strength for isotopes below the mass 
A = 96, which is not surprising because our model space 
allows excitations out of the pi/2 shell. On the other 
hand, it overestimates the strength for some isotopes with 
A = 97, which is related to the partial occupation of 
the 37/2 shell by neutrons reducing the total strength as 
compared to the non-interacting picture. 

Since our calculated total strength is reasonably close 
to the single-particle estimates, we could use the exper- 
imental hindrance factor quoted relative to the single- 
particle estimate: h'^^P — GTgp/GTexp- Unfortunately, 
in this region, experimental information on h'^^^ is lim- 
ited. Some of the calculated nuclei naturally decay from 
the ground state by /?" -decay instead of electron cap- 
ture, while in other nuclei, the Q-window contains only a 
small fraction of the total strength. Thus the total GT+ 
strength could be obtained only by (n,p) or similar one- 
particle exchange reactions. An additional uncertainty in 
deriving the hindrance factor, even for nuclei where the 
Q-window is large, comes from the fact that 7-ray spec- 
troscopy misses a significant fraction of the Gamow- Teller 
decay strength due to sensitivity limits of detectors, a low 



population of nuclear levels close to the Q-limit as well 
as weak intensity of their decays [10, 21]. This limitation 
can be largely overcome combining a high-resolution 7- 
ray detector with total absorption spectrometry (TAS), 
as was done in a number of recent experiments on nuclei 
in the ^°"Sn region. For example, a study of ^^Ag de- 
cay [11] showed that only 2/3 of the total Gamow- Teller 
strength is obtained by high-resolution 7-ray spectrome- 
try, while the same number for ^''^In was about 1/8 [21]. 

One nucleus, ^^Ag, has almost 98% of the total GT 
strength inside the Qec window. We can use this 
nucleus to estimate the experimental hindrance factor. 
Hu et al [11] reported Y.B{GT) = 3.00(40) based on 
TAS measurements, which leads to the hindrance factor 
li^xp _ 4,24to'5o- Another nucleus where this window is 
large, and there is a TAS measurement available is ^*Ag 
[22]. Hu et al reported the total strength in ^*Ag to 
be 2.7(4), giving the hindrance factor h^^'P = 4.27l:J5;^5, 
since the calculated Bqj^^. is 11.53 (this is 92% of the to- 
tal Gamow- Teller strength inside the Qec ~ 8.24 MeV 
window). Thus the hindrance in two Ag isotopes, ^^Ag 
and ^^Ag, are of the order of 4.25, and we adopt this 
value for the total hindrance factor h. We did not con- 
sider heavier nuclei for the hindrance estimate, because 
they are further away from our region of interest, and the 
possible Z-dependence of this factor is not clear [11]. 

In this region, the only available total Gamow- Teller 
strength measured using {n,p) reaction is for ^°Zr. Ray- 
wood et al [23] deduced a value of 1.0 ± 0.3 for the total 
strength. Our calculated total strength for this isotope 
is S{GT+) = 0.34. We should note, however, that in our 
restricted model space there is only one 1+ state in ^"Y. 
These values can also be compared to recently reported 
measurements of 3.0 ± 1.9 for the total strength by Sakai 
and Yako [24]. 

We turn now to strength distributions. Our calculated 
Gamow- Teller distributions in the decay of nuclear sys- 
tems with a few valence protons, like Zr or Mo isotopes, 
have the strength concentrated in a narrow energy range 
(less than 0.5 MeV) or sometimes in only one transition. 
The strength in systems with ZJ > 4 (Tc and above) is 
distributed over the energy range of about 4 MeV. The 
Nb isotope chain is intermediate in this respect, because 
in the lowest configuration Nb has only one valence pro- 
ton in the 59/2 shell, while decays to Zr isotopes are dis- 
tributed over several states. We show these systematics 
using a few examples below. 

In Figs. 3-5, we compare the calculated Gamow- Teller 
distributions with available data collected from several 
sources. All measurements were done using 7-ray spec- 
troscopy. Data for A^ = 50 and 51 isotopes were obtained 
from Refs. [10, 25, 26] (see also references there in). The 
GT distribution in ^^Ag [11] was obtained with a TAS 
measurement. We note that in comparisons to experi- 
mental data, we do not include the sensitivity limits of 
experimental detectors. This sensitivity artificially cuts 
off Gamow- Teller strength near the Q-window so that 
calculated states are often not observed even 2 MeV be- 
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FIG. 3: (Color online) Gamow- Teller strength distributions 
from the ground states of A'^ = 50 isotones as a function of 
the excitation energy in the daughter nucleus. Two experi- 
mental sets are plotted for the decay of ^^Rh (the values are 
taken from [10, 25]). The gray-shaded area is the calculated 
strength, while the histogram represents experiment. The ex- 
perimental data suffer from detector sensitivity limits near the 
Q-limit (indicated by an arrow), with an exception of ^^Ag 
where the data were obtained using TAS [11]. We show the 
strength in ^'^Ag scaled by a factor of 2. The bin size is 0.2 
MeV. 
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FIG. 4: (Color online) Gamow- Teller strength distributions 
from the ground states of two A^ = 51 isotones: ^^Ru and 
^'^Pd (values are taken from [26]). See also the caption of Fig. 
3. 
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FIG. 5: (Color online) Gamow- Teller strength distributions 
from the ground state and first excited states of ^®Rh (values 
are taken from [26]). See also caption of Fig. 3. 



low the Q-window (see, e.g., [10, 11, 21]). For this rea- 
son, coniparisons to experiment are somewhat difficult, 
and one should focus on unambigous regions of low- lying 
strength. In Table I we listed fractions of the calculated 
Gamow- Teller strength that lie inside the Q-window. A 
comparison to the calculation by Brown and Rykaczewski 
[10] reveals some interaction dependence in the values. 
For example, they estimated Jec — 29% and 99% in the 
decays of ^''Ru and ^^Pd, while our estimates are only 
19% and 72%, respectively. Johnstone [25, 26], following 
a different approach, estimated significantly higher frac- 
tions of the strength inside the Q-window for most cases 
except ^^Ru. 



TABLE I: Fraction of the calculated Gamow- Teller strength 
inside Qec window, /ec = Bqj^^/GT+ 



Reaction 


GT+ 


Bqec 


/bc, % 


93Tc(/3+)93Mo 


6.12 


0.13 


2 


S4Tc(/3+)34Mo 


5.74 


0.13 


2 


s*Ru(/3+)94tc 


7.89 


1.53 


19 


95Tc(/3+)95Mo 


5.43 


0.08 


1 


^'Ru(/3+)9STc 


7.49 


0.59 


8 


35Rh(/3+)95Ru 


9.41 


6.35 


67 


96Tc(/3+)9«Mo 


5.36 


0.04 


1 


s«Rh(/3+)*Ru 


9.05 


5.14 


57 


96pd(/3+)*Rh 


11.18 


8.08 


72 


9^Rh(/3+)9^Ru 


8.53 


2.60 


30 


s^Pd(/3+)97Rh 


10.90 


7.38 


68 


^'HiP+rVA 


12.71 


12.49 


98 


9«Ag(/3+)3»Pd 


12.47 


11.53 


92 



From Figs. 3-5 we note that the calculated strength 
distributions follow the trend observed in experiments. 
Most odd-Z A^ = 50 isotones and A^ = 51 isotones have 
little strength at low excitation energies [E^ < 2-3 MeV) , 
with the strength distributed among many states at a 
higher excitation energy, some of which are above the Q- 
window. The strength in even-even nuclei, represented 
here by even-Z A^ = 50 isotones, is concentrated in a few 
states. 

The Gamow- Teller distribution in ^^Ag shown in Fig. 
3 is converged (around 60 iterations per J/ was required); 
thus the calculated shape is as good as it can be for the 
interaction. The ccntroid of experimental Gamow- Teller 
strength distribution in ^^Ag is lower than the calculation 
predicts: El:'J^^^ = 4.3 MeV versus E^^f^^^ = 4.7 MeV. 
This is one of the indicators that the interaction may 
require some fine-tuning. 

We already showed and discussed parts of the calcu- 
lated Gamow- Teller distributions in Figs. 3-5. Having 
in mind upcoming experiments [9] on Mo isotopes, we 
show the calculated distributions for the decays of Mo 
isotopes with masses A — 93-97 in Fig. 6. The main con- 
tributions to the total strength are located within 1 MeV 
energy range around the centroid. The decays of Tc iso- 
topes (see Fig. 7) have the strength distributed within 4 
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FIG. 6; Calculated Gamow- Teller distributions in Mo iso- 
topes with masses A = 93-97 (A^ — 50-55). The arrow indi- 
cates that the value inside the bin (the value is shown by a 
number beside it) is greater than the upper limit of the j/-axis. 
The bin size is 0.2 MeV. 
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FIG. 7: Calculated Gamow- Teller distributions in Tc isotopes 
with masses A — 94-97. See also caption to Fig. 6. 



MeV range. These two isotope chains display the differ- 
ence in the decays of even-even and odd-odd nuclei that 
we discussed above. 

We also show the distributions for Zr isotopes, see Fig. 
8. There we observed the migration of the strength from 
higher energies to lower energies as the neutron number 
increases. At ^^Zr, where the lowest configuration is a 
completely occupied proton pi/2 shell and a similar situ- 
ation occurs in the neutron d^/2 shell, the entire strength 
is peaked in one transition. The trend remains also in 
97Zr. 

We turn now to a discussion of the calculated total 



FIG. 8: Calculated Gamow- Teller distributions in Nb isotopes 
with masses A — 92-96. See also caption to Fig. 6. 



Gamow- Teller strength and the centroids. The isotopic 
dependence of the total strength is smooth with the 
strength decreasing together with the increasing num- 
ber of neutrons and/or the decreasing number of valence 
protons in the 179/2 shell. Assuming no mass depen- 
dence, an approximate formula can be derived: GT^ = 
0.086(Z^ - 1.5)(20 - 7V„). (The factor (20 - N^) is due 
to the relative unimportance of the /in /2 shell because 
of its negative parity.) This form is somewhat similar to 
the dependence Zy {20 — Ny)/A observed in the pf shell 
nuclei (see e.g. [27]). The difference may be related to 
the active j shells. In the pf shell nuclei, discussed in 
Ref. [27], protons predominantly occupy the f'j/2 shell; 
thus its occupation is proportional to the number of va- 
lence protons. While in our model space, the occupancy 
of the gg/2 proton shell increases due to excitations out 
of the pi/2 shell via configuration mixing. This increase 
is greater for isotopes closer to the core (around 0.6), and 
is 0.2 for A = %7 nuclei with Zy > 2. The formula's y^ 
per degree of freedom is 0.05. 

Another systematic relates to the centroids of the GT^ 
distribution. If plotted with respect to the lowest posi- 
tive parity state of the daughter nucleus, the centroids 
of Gamow- Teller distributions show a characteristic odd- 
even staggering, see Fig. 9. They are low in even-even 
nuclei, high in odd-odd nuclei, and average in odd-^ nu- 
clei. A similar trend was observed in the vaid-pf shell 
nuclei [28]. 

Langanke and Martmez-Pinedo [5] interpreted this 
odd-even staggering as a result of pairing energy con- 
tributions to the mass splitting between the parent and 
daughter nuclei (see also [29]). The pairing structure 
goes away if the centroids are measured with respect to 
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FIG. 9: Excitation energy of the centroid in the calculated 

nuclei. For Zr-^Y decay, the excitation energy was calculated 
from the lowest positive parity state in Y. 
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FIG. 10: (Color online) Energy of the GT+ centroid with 
respect to the calculated ground state energy of the parent 
nucleus as a function of (A'^ — Z)/A. See also text. 



the parent nucleus. We plotted the centroid energies cal- 
culated in this way in Fig. 10. We also included the 
Coulomb energy difference, calculated using a formula 
[30] Ec = 0.72(Z7Ai/3)(l - 1.69M2/3), but ignored the 
proton-neutron mass difference and the splitting between 
the proton and neutron single-particle orbits which would 
be present if the lowest single-particle energies would be 
taken with respect to the core nucleus, ®*Sr. The fig- 



ure shows that centroid energies indeed lose information 
about the pairing structure. It is interesting to note that 
there seems to be a cross-over behavior, which we high- 
lighted by connecting the transitions corresponding to 
the decays of Zr isotopes. These centroid energies follow 
a linear dependence as well, but the inclination is differ- 
ent from that in other nuclei. This behavior is probably 
related to the fact that GT+ strength in Zr is due to 
proton excitations out of the Pi/2 shell. 



IV. DISTRIBUTED-MEMORY SHELL-MODEL 
CODE 

Our calculations were performed using a new parallel 
shell-model code ORPAH (Oak Ridge PArallel shell model 
code), which is still under development. The basic ideas 
are similar to those employed in the serial m-scheme com- 
puter code ANTOINE [31, 32]. However, there are differ- 
ences in the approach, since the code was developed tar- 
geting the distributed-memory computational paradigm. 
While the distributed-memory approach sets no limits 
on the available memory or the number of processors in- 
volved in the computation, a natural limitation occurs 
due to the need to communicate data from one processor 
to another, a process which for collective operations scale 
as iVp, where Np is the number of processors. However, 
even in cases when the communication becomes unfavor- 
able, there is still a possible trade-off because of a greater 
amount of available memory. 

The most time-consuming part of the shell-model prob- 
lem is the operation of the Hamiltonian on a vector to 
produce a new vector which occurs during the Lanczos 
procedure. This affects load-balancing since each pro- 
cessor needs a similar workload for an effective use of 
computational time. We consider parallelization at two 
levels: the vector amplitudes are distributed among the 
processors, and each processor produces a portion of the 
final vector in a time-balanced way. We then redistribute 
the final partial vectors to the appropriate positions after 
each iteration. 

Similar to ANTOINE, the code numerically builds 
"blocks" of identical-particle Slater determinants having 
the same quantum numbers and sets up tables allowing 
construction of the elements of the Hamiltonian matrix 
[32]. Differences arise from the parallel implementation. 
If the dimension of the model space is £), then the part 
of amplitudes which reside on a particular processor has 
size D/Np. For a sufficiently large number of proces- 
sors, D/Np can get smaller than the size of the largest 
block. This would place the limit for the maximum rea- 
sonable number of processors. Operations involving this 
block are also the most time-consuming. Wc decided to 
split the blocks in order to have smaller pieces of tasks, 
which could be distributed among the processors more 
efficiently. During the operation of the Hamiltonian act- 
ing on a vector, some of the amplitudes are prefetched 
and others are requested during the calculation. There is 



no predefined communication pattern, since some ampli- 
tudes can be delivered with some delay while the process 
would still be able to employ the ones already present in 
the memory. To enable this disconnection of computa- 
tion and communication, each processor consists of two 
threads responsible for those two tasks. Those threads 
communicate via shared variables. The need to deliver 
amplitudes creates a communication overhead on top of 
the time needed to produce the final Lanczos vector. The 
latter operation is well balanced (i.e., is inversely propor- 
tional to the number of processors), while the overhead 
depends on the number of processors involved in the cal- 
culation. The final Lanczos vector is reorthogonalized to 
previous vectors. Due to orthogonality of the basis, each 
processor can produce the partial sum of a scalar prod- 
uct, and a global communication is needed only to obtain 
the total sum. 

There is no coded-in restriction on the number of pro- 
cessors, with the exception that the minimum number of 
processors is two, because of the manager-worker algo- 
rithm employed in the Hamiltonian table set-up proce- 
dure. The current version of the code can calculate eigen- 
values and eigenvectors of the Hamiltonian, the total an- 
gular momentum and isospin, as well as GT properties. 
Some computations were done on a 2-4 CPU computer; 
others were done at the NERSC computer Seaborg using 
up to 80 processors. The largest problem that we tried to 
solve was the ground state energy of ^^Fe (Z) = 110 x 10^) 
on 48 processors, but the current set-up did not allow 
us to reach such dimensions in the region of our study. 
In addition, a further improvement in the performance 
is required before the code assessment is done, though 
the distributed-memory computation is a venue to solve 
larger interacting shell-model problems. 



SUMMARY 



values slightly differ from a naive single-particle shell- 
model picture. Calculated strength distributions appear 
to reasonably recover experimental distributions in re- 
gions that are unaffected by detector sensitivity limits. 
From TAS data on ^^Ag, we were able to obtain an esti- 
mate of the phenomenological quenching factor relative 
to single-particle estimates. Furthermore, our GT distri- 
bution for ^^Ag reproduces the measured data, see Fig. 
3. By analyzing the centroids of the Gamow- Teller dis- 
tributions, we found that the odd-even staggering be- 
havior disappears if the centroids are measured from the 
parent ground state, as was suggested by Langanke and 
Martincz-Pincdo [5] . We also observed that the centroids 
measured in this way have a quasi-linear dependency on 
the parameterization (N — Z)/A, with a different inclina- 
tion for nuclei where no protons are present in the gg/2 
shell in the non-interacting picture. Finally, we made 
several predictions of the strength distributions for mea- 
surements that may soon be available in mass A — 92-97 
nuclei. 

The description of low-lying Gamow- Teller strength 
distributions in a large range of nuclei (from roughly 
mass A ~ 50 to mass A ~ 150) is one important in- 
gredient in understanding type-II supernova explosions, 
since electrons get captured by nuclei through these lev- 
els. Of course, this is not the whole story since first- 
and second-forbidden transitions (typically difficult to 
access in the laboratory) also play an important role in 
the cross section and rate calculations relevant for super- 
nova. For low-energy capture, the Gamow- Teller tran- 
sitions will dominate. They also dominate the neutrino- 
nucleus scattering processes that may occur at later times 
in the supernova event. While theory can produce many 
things, measurements are necessary to validate the esti- 
mates. Such is the case in the nuclear region discussed in 
this paper. We eagerly await experimental comparisons 
to our calculations. 



We calculated nuclei above **Sr having masses A = 90- 
97 using a realistic effective interaction derived from the 
CD-Bonn potential. The agreement between the calcu- 
lated and measured spectra is satisfactory. Improvements 
to the interaction through fine tuning of the matrix el- 
ements could be useful to obtain the finer spectroscopic 
details, including the level ordering or the placement of 
negative-parity states in several nuclei, but this is beyond 
the scope of this exploratory work. 

We also calculated the total Gamow- Teller strength 
and strength distributions for the decay in the electron 
capture direction. We found that the total strength 
follows the single-particle estimate based on the TTgQ/2 
and i'gr/2 occupation numbers obtained from the ground 
state wave functions of the parent nucleus, although the 
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